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1> ' Abstract 

00 . 

We present results of an investigation into the nature of instantons 

in 4-dimensional pure gauge lattice SU(2) obtained from configurations 

which have been cooled using an under-relaxed cooling algorithm. We 

iy~j . discuss ways of calibrating the cooling and the effects of different degrees 

of cooling, and compare our data for the shapes, sizes and locations of 

instantons with continuum results. In this paper we extend the ideas and 

techniques developed by us for use in 0(3), and compare the results with 

Q_i| those obtained by other groups. 
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1 Introduction 

Much work has been done in simulating the SU(2) vacuum, with the aim of 
learning more about the vacuum structure of QCD, and a lot of interest has 
been shown in the existence of instantons in this vacuum. Early studies into 
the cooled SU(2) vacuum were carried out in |L|, where, working on a lattice 
with 'zero-derivative' boundary conditions rather than the more usual periodic 
boundaries, Teper investigated the role of instantons by locally minimizing the 
action (ie. cooling) and concluded that the lattice SU(2) theory does indeed 
contain instantons, but that the large, physically important ones will be absent 
in studies with (3 < 2.7. At the same time, in [F2J, Ilgenfritz et al, working on 6 4 
lattices at f3 = 2.1, 2.2, found (multi-)instanton solutions from successive coolings 
of the (equilibrated) gauge fields: 

U x ^ — > U' X fl = cL x ^ (1) 

where c is a normalisation constant to ensure „ G SU (2) and E is the sum of 
six staples (this will be equivalent to taking a = in eq. |2"T|). They found that 
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the action density is concentrated in distinct, separated parts of the lattice - 
consistent with the picture of a dilute instanton gas — and that these lattice 
instantons are meta-stable under cooling of the gauge fields, being shrunk in size 
and finally annihilated by prolonged cooling. They also studied the eigenvalue 
spectra of the fermion zero modes, and concluded that £77(2) on the lattice does 
indeed have an instanton vacuum. 

In ||, Polikarpov and Veselov, working on 8 4 -12 4 lattices at /3 — 2.3,2.35 
carried out a detailed study of the cooling process and identified five distinct 
stages in the history of a cooled configuration: (i) freezing of quantum fluctuations 
and the arising of quasi-classical configurations of instanton-anti-instanton type; 

(ii) smoothing of quantum fluctuations around instantons and anti-instantons; 

(iii) annihilation of instantons and anti-instantons; (iv) a long stage involving 
the existence of classical solutions; (v) a short stage involving the disappearance 
of classical solutions. 

The classification of these stages is, in their own words, somewhat arbitrary, 
with an absence of definite borders between the stages. They found that there 
were (at least) two characteristic points on the cooling curve, however: around 
the middle of stage (iv) they found a zero of A 2 S(N) (corresponding to the second 
derivative of S with respect to cooling sweep N), and another somewhere on the 
border between stages (ii) and (iii) at which A 2 S(N) « 0. They calibrated their 
cooling process by stopping at the former point; they denoted a configuration 
in this state as an 'instanton vacuum', and identified the latter point as the 
beginning of the instanton-anti-instanton annihilation, denoted as the 'instanton- 
anti-instanton vacuum'. 

This paper is organised as follows: in section we outline the formalism of 
continuum SU(2), and briefly discuss the instanton solution found by Belavin 
et al in @]; in section || we consider how to construct an SU(2) instanton on 
a periodic lattice; in section f| we discuss calibrating the cooling we use, with 



a first calibration attempt given in subsection fO| and a better method given 



in subsection O; in section |5| we detail the means of locating and sizing the 



instantons, compare their shapes with the form of isolated, continuum objects 
and calculate a size distribution for them; in section ^| we consider inter-instanton 
separations and compare our results with refs. 0. Our conclusions are given 
in section |7|. The figures referred to in the text are collected together at the end 
of the paper. Preliminary results from this work were presented at Lattice 94 [7]. 



2 Instantons in SU{2) 

As in 2-dimensional 0(3), 4-dimensional SU(2) has instanton solutions^ due, in 
this case, to the homotopy class of mappings from S 3 — > S3, ie from the points 

1 This section follows the discussion of Yang-Mills instantons given in M. 
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on a four-dimensional sphere parameterised by three angles, to the elements of 
SU(2), which is also determined by three parameters. In other words, the SU(2) 
group manifold is topologically equivalent to S3: any SU(2) group element can 
be written in terms of the Pauli matrices as 

U = m 4 1 + iu ■ er ; u\ + u 2 = 1 , er = (cti, cr 2 , 0-3) (2) 

with 1 the 2x2 unit matrix, and u = (u^, u) clearly defining a point on S3, the 
sphere in 4-dimensional Euclidean space. Mappings in this case are characterised 
by the topological charge Q T defined in terms of the divergence of a (gauge- 
dependent) current: 

Q T = -77n, I d * x tr F ^ ( 3 ) 

= 2tr (F^) (4) 

2 

K a = 4£ a(hS tT {ApdyAg + -ApA^As) (5) 

with the gauge fields and field strength tensor scaled by i/g and defined by 

= d„A v - d v A„ + [An, A v \ (6) 

and the dual of F given by F^ v = \spuX P Fx P - Under a gauge transformation 
T(Xfj) we have 

A^T^A^T + T^dpT . (7) 

To find the instanton solutions we look for finite-action solutions to classical 
Euclidean Yang-Mills theory. In order to have the action 

1 



S = — J d*x tr F^ (8) 

finite we require F fJ-l/ (x) '^t^ 00 which we can take to mean A^x) '^U* 00 0. 
However, in practice this is too restrictive, and eq. |7| means we only need 

Ap(x) -> T^dpT as \x\ -> 00 (9) 

which is obtained from A^ = by a gauge transformation. To find the fields 
satisfying this boundary condition we make use of the following positivity relation 



tr I d x (Fp V ± F^v) — • (10) 

As 



(F^u ± Fp V f = 2(F flu F /lu ± F^F^) (11) 



we have 



tr / d*x F^F^ > 



tr I d x FpyF^y 



16tt 2 Q (12) 
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using the definition for Q given in eq. |3]. From this we have 



S > KQ (13) 
T 

whence we see that the action for a single instanton is Si = 87c 2 /g 2 (cf. the result 
for 0(3): Si = An/g 2 ). Clearly, then, the action is minimised for the (anti-)self 
dual solutions F^ u = ±F^ U . 

In looking for the instanton solution, Belavin et al M considered first 0(4) 
gauge theory, isomorphic to SU(2) x SU(2) with one SU(2) identified with the 
internal symmetry and the other with the three-sphere at \x\ — > oo. The resulting 
gauge transformation has the form 

T(x) = -^-(^ + ix ■ <x) with \x\ 2 = x\ + x 2 (14) 



\x\ 



giving rise to a gauge field 



x 2 



A,(x) = - T —- 2 T^(x)d fl T(x) (15) 
p l + x 2 

where p is the instanton size, and x = its centre. For x 2 3> p 2 , — > T^d^T as 
required by eq. ||. We can write the space and time components of explicitly 

as 

A 4 (x) = -^—^ (16) 
x z + pr 

x z + pr 

2.1 Instanton— ant i- instant on interactions 

Thus far we have only considered the single-instanton solution; multi-instanton 
configurations can be generated by simply multiplying instanton fields with their 
centres at different locations |J , and calculations involving mixed instanton-anti- 
instanton configurations are non-trivial. In ||], Forster considered an instanton- 
anti-instanton pair at large separation, taking an instanton at x = 0, given by 



eq. 15, and an anti-instanton at x — C, given by: 



= x 2 + (x°- C) 2 T{X ~ °) d ^ x - C ) ( 18 ) 

with A the size parameter for the anti-instanton. By performing a conformal 
mapping, the point x = C is sent to infinity and the local behaviour at x = C 



4 



implies the interaction takes the form of a kink and an anti-kink separated by a 
distance — log(pA) giving an interaction term: 

Sint = -48^ p2>? (19) 

with g 2 the coupling occuring in eq. |8| above, implying an attraction between 
unlike objects. In ||, Palmer and Pinsky, investigating this interaction and the 
distinction between dilute and dense phases of the instanton gas, considered the 
case of an instanton-anti-instanton pair separated by the size of one member, 
with the size of the other free. They found that there was a minimum finite 
separation: 

c = p = 2A (20) 

at which point the pair melts into a pure gauge with zero action. This will be 
relevant when we investigate the separations of objects in section |6] below. 

3 SU(2) instantons on the lattice 



We seek to extend the successful calibrated cooling method [10, 11 1 from 0(3) to 



SU(2). As before, we use an under-relaxed cool, in this case given by 

U'^^E^ + aU^ (21) 
where is the sum of the six 'staples' around the link U X)ll , and U' is sub- 



sequently normalised to lie in SU(2). Following the results of [|10], |J we take 
a = 2. 

In order to calibrate the cooling, we want to generate a gauge configuration by 
using a suitable lattice version of eq. |T|, as we intend to calibrate the cooling by 
examining its effects on a particular known configuration. However, this formula 
is explicitly incompatible with periodic boundary conditions^, as we require the 

fields to vanish at large distances from the instanton centre^]. In order to 
achieve this we must construct the instanton as follows [|T^]: first construct the 
A^(x) in R A as above. Next construct the link variables U XifM from these continuum 
fields. As stated above, we require an instanton solution to be self-dual, and for 
the lattice solution this relation is only approximate. If we take as the link fields 

U Kili = exp(i[A^x) + A^{x + n)\/2) (22) 

constructing these on an infinite lattice, then the deviation from self-duality is 
0(p~ 2 ) |12|]. Next we perform a singular gauge transformation on the lattice 

1 We are greatly indebted to Mike Teper for discussion on this subject. 

2 The concept of "large distance" on a periodic lattice may seem a strange one. What we 
need, for a moderately sized instanton sited at the centre of the lattice, is that the fields 
approach zero at the lattice boundary. 
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using eq. |14], with the singularity 'off-site' both here and in the construction of 
the A Xt/Ji fields: 

U x ,fj, = T x U X) pT x+ ^ (23) 

(where T x denotes the gauge transformation on the site x of the infinite lattice) 
and finally we impose periodicity on the solution by enforcing 

U^Xi = -L/2) = U^Xi = +L/2) (24) 

for each Xj and each /i for a L 4 lattice. 

In constructing this solution we must also take into account the fact that 
eq. |1| gives the instanton in terms of the A^(x) fields, while on the lattice we 
need to use the U x<tl link fields, and so need to rewrite the instanton solution in 
terms of these link variables. This is straightforward as 

exp(i6*a; • cr) = cos8 + i(x ■ cr) sin 9 (25) 

and we can rewrite eq. [T5]m this form. The individual components are, explicitly: 

L4 4 (x) = (cr ■ (x 1 ,x 2 ,x 3 )) (26) 

x z + p z 

L4i(x) = 1 (cr ■ (x 4 ,x 3 , -x 2 )) (27) 
x z + p z 

iA 2 {x) = (sigma - (-x 3 , x x )) (28) 

x z + p z 

iA 3 (x) = (sigma • (x 2 , -x h x 4 )) . (29) 

x z + p z 

As in 0(3), this field configuration is metastable (with the Wilson gauge action 
we are using) and will be annihilated by prolonged cooling. (See the discussion 
in (1^, [ll[). For instantons where p is not small compared to L there will be 



significant discontinuities across the lattice boundary, which can be removed by 
a moderate amount of cooling. 



4 Calibrating the under-relaxed cooling. 

We performed our simulations on a CRAY Y-MP, using 16 4 and 24 4 lattices at 
f3 = 2.4, 2.5 respectively. These correspond to lattice spacings of approximately 
0.12fm and 0.08fm ||13| . Each of our configurations is generated from 500 up- 
dates with 1 heatbath update followed by 3 over-relaxation updates, after several 
thousand thermalisation sweeps. In fig. [I] we present the cooling history of two 
example SU(2) gauge configurations, calculated on the smaller lattice. 
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4.1 A first attempt at calibration 

For 0(3), we concluded JTTJ] that an optimum cooling would be that required to 
annihilate an instanton of size p = 0(2a) where a is the lattice spacing. In table | 
there is shown the number of cooling sweeps at a = 2 (which we had decided was 
a good compromise between cooling too severely and making too many cooling 
sweeps) needed to remove SU(2) instantons of size p = a, 3a/2,2a,4a generated 
from the prescription above. It transpires that these are much more stable under 
cooling than their 0(3) counterparts. 



p a 3a/2 


2a 


4a 


iV sweeps 495 503 


510 


536 



Table 1: The number of sweeps N sweeps at a = 2 required to annihilate (to reduce 
S/Si to < 0.1) a single SU(2) instanton of size p. Calculated on a 16 4 lattice at 
(3 = 2.4, with the instanton positioned off-site in the centre of the lattice. 

This much greater stability presents us with something of a dilemma: were 
we to naively proceed and take the value for p = 2a in table [I] as a calibrator, 
then the gauge configurations we would have after cooling would be well past the 
first three stages defined by Polikarpov and Veselov above — which is where we 
wish to examine the vacuum, as we wish to study its structure in the presence 
of (anti-)instantons - - and in some cases will have been cooled so far as to 
entirely remove all vacuum structure. There are two possible candidates from 
the cooling curve of annihilation of a single instanton for the level of cooling to 
use. The first is the cooling required to remove the undesired boundary effects, 
while the second is the amount required to then annihilate the instanton, once 
the boundary discontinuities have been eliminated. However, neither of these 
methods of calibration is really satisfactory — the former is too unphysical a 
process, as it removes what is essentially a lattice artefact that has nothing to do 
with the underlying phyics, and given the size of the lattice and the value of (3 
we use, the latter seems excessive when compared to the results presented in [ |TT|] . 
In light of this we propose a different means of calibration. 

4.2 A second calibration method 

As we wish to model the vacuum as quantum fluctuations around classical config- 
urations, we first generate a classical configuration: the single-instanton solution 
as presented above. We then cool this by enough sweeps to remove any boundary 
discontinuities. Next we perform a few heatbath updates at the (3 value under 
study to introduce some quantum fluctuations and then cool it until we regain 
the original configuration, ie until S — Si ± 10%. 
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Quite how many heatbath sweeps to perform is open to question: we seek to 
reach a state in some respects similar to the uncooled configurations, but not too 
far away from the original instanton — we do not want to introduce any more 
structure than was already present. We found that, at (3 = 2.4 on a 16 4 lattice, 10 
heatbath sweeps was optimum for introducing the required quantum fluctuations 
without overly changing the underlying structure, whereas for (3 = 2.5 on a 24 4 
lattice, 10 was too few and, at this higher (3 value, we found 20 were required. As 
a check of this we found 20 too many for the lower f3 on the larger lattice, with 
new structure occuring in more than half of the sample configurations we looked 
at; by comparison, for (3 = 2.5, 20 sweeps introduced new structure^ in less than 
5% of the ensemble. 

As a further check of the validity of the number of thermalisation sweeps used, 
we looked at the value of (trf/ n ) on the configuration. We found that with 10 
heating sweeps on the (3 = 2.4, 16 4 lattice this was consistent with established 
results, and this was also true for 20 sweeps on the (3 = 2.5, 24 4 configurations, 
while half the number of sweeps in each case gave a value that was too high. These 
two checks — not introducing more structure, and heating until the correct value 
of (trU n ) is attained at a given (3 value — form respectively effective upper and 
lower bounds on the number of heatbath sweeps to be performed. 

The calibration results are given in table ^| for initial configurations each 
containing a single instanton of size p = a, 2a, 3a, 4a. 

In [[ll|] we decided that a good check on whether the cooling had taken the 
configuration to the same physical state was whether or not there was a good 
agreement across (3 of the value of (SyV^hysSV), and found that the criteria we 
used for 0(3) led to such an agreement. For our SU(2) data however, the agree- 
ment appears to be less good, with (S/V p hy S Sj) = 0.000185(2) at (3 — 2.4 and 
(S/V phys Si) = 0.000249(3) at (3 = 2.5. While not equal, these two numbers 
seemed sufficiently close for us to accept this calibration method, although no 
doubt these values could be brought closer by a certain amount of 'fine-tuning'. 
On the basis of these results we chose 72 sweeps to be the number required at 
(3 = 2.4 with a = 2, and 154 sweeps at (3 = 2.5, and accepted this method as a 
calibration technique for under- relaxed cooling in 4d SU(2). 



5 Instanton size, shape and distributions 



In ||TT| we found we could obtain good results for the sizes and locations of 
instantons using S(x), although in principle Q(x) contains more information. 
This was because in investigating the size distribution we were not concerned with 



■""We ascertained whether the configuration had been moved away from the single-instanton 
solution by following the cooling history. Where new structure had been introduced by the 
heatbath sweeps, 'steps' appeared in the history, in much the same way as they appear in the 
solid line in fig. [j]. 



s 



(5 L iV heat p N cool (trUu) 



2.4 16 10 a 


71(8) 


0.6344(6) 


2a 


72(7) 


0.6346(5) 


3a 


71(7) 


0.6347(7) 


4a 


73(5) 


0.6345(5) 


2.5 24 20 a 


152(12) 


0.6529(3) 


2a 


155(8) 


0.6530(3) 


3a 


154(11) 


0.6530(3) 


4a 


154(13) 


0.6530(3) 



Table 2: The number of sweeps iV coo i at a = 2 required to restore a config- 
uration with a single lattice instanton of size p after iV heat heatbath updates, 
and the value of the average plaquette after the heatbath sweeps. The figures in 
parentheses are one standard deviation, the data are obtained from 50 configu- 
rations in each case. The accepted values of (trZ7 a ) are 0.63058(2) [j3 = 2.4] and 
0.65235(l)[/3 = 2.5] 0. 



the sign of the object. When we came to look at the inter-instanton separations, 
we still used S(x) for calculating the locations and obtained information on the 
sign from Q(x). 



For SU(2), we used the 'plaquette' method (T!J to calculate Q(x) on our 
cooled configurations: 

= —r^ (tr (n Xil2 n x , 34 ) + tr (IV3IV2) + tr (n^IL^)) (30) 

where H X) ij is the orientated path combination at x in the ij plane, consisting of 
the sum of the four clover leaves of (U a — U^)/2. 

We choose the most local and symmetric definition for the action density at a 
site. We used the average of contributions from all plaquettes with one corner on 
the site x. This results in an action density that matches Q(x) very well (and in 
the continuum S(x) oc Q(x) around an isolated instanton), but is not as smooth 
as Q(x), as can be seen from figs. |2| and [| Since S(x) not as smooth as Q(x), on 
occasion the instanton location algorithm will record more structure in a given 
area of the lattice than is actually there and will ignore valid structure elsewhere 
when applied to S(x). For this reason we conducted our instanton search at 
/3 = 2.4 using both S(x) and Q(x). 

Here we look at connected 4-dimensional regions of the lattice. We use the 



same method as in |TT| by looking for local maxima in the density, and in this 
case use as a size parameter, which we denote pL, the 4th root of the hypervolume 
of the connected region of sites for which the density is not less than half that 



9 



at the peak. In order to relate this pl back to the size parameter in eq. [15], we 
define an effective radius (see ref. |16| for a similar prescription) p eff : 

2 2 /oi\ 

P '« = ^W^V) PL - (31) 

This is the parameter plotted in figs. [§-[7| and quoted in table |3|. Eq. [13] gives a 
lower bound on the action, and so S/Sj can be interpreted as an upper bound 
on the number of objects in a given configuration, and we used this to determine 
how many we looked for in each configuration.. 

To see how closely the individual objects resemble isolated continuum instan- 
tons, we looked at the relation between the height at the peak and the calculated 
size, as these have the continuum relations 

'S'max = — o~ a (32) 

9 P 



O --=^-9 -— — (W) 

^fmax 1e 2 max 2 4 V / 



and 

2<7 2 c _6_ 

167T 2 ^ 7T 2 P 4 

for an isolated instanton. Our results are presented in figs. f| and |[ All distances 
are in units of the lightest glueball mass, £ = l/m + with m +(2.4) = 0.94(8) and 
m +(2.5) = 0.63(8), taking the values for the lattice spacing a at these values of 



(3 from |T3| and the mass of the + glueball at (3 = 2.4 from [17]. The data was 
obtained from 261 configs at (3 = 2.4 and 250 at (3 = 2.5. 

As in |ll|] we calculated the distribution of instanton sizes. The data for 
dN/Vdp are presented in figs. |6|-[7[ All sizes are again in units of the correlation 
length, £. (It should be noted that the plot we presented earlier in [0 used p^, 
scaled by a different factor from that above, rather than p eff and so contained a 
peak in a different position.) 

In agreement with 0, [3], |18| we find a distribution that grows with instanton 
size (unlike that presented for 0(3) in |nj), although our distribution lacks the 
sharp large-size cutoff they find; we attribute this to our means of determining 
the size from connected hyper-volumes whereas in [§, || they simply measured 



the size along the lattice axes. In ||18|| , Chu et al measured the size by using the 
correlation function 

f(x) = J2Q(y)Q(x + y) (34) 

y 

and comparing this with a convolution of the analytic expression: 

Qp( x ) = -yi (-ax-*) ( 35 ) 

7r 2 /r \x 2 + p l j 

which they fitted to the lattice data for f(x), obtained by cooling for 25 and 
50 sweeps, using a single value of p they took to be an average size. Their 
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distribution was obtained, like ours, by analysing connected regions of the lattice, 
and extracting the size from the size of the region, their criterion being that two 
adjacent points belong to the same cluster if the product of Q(x) at these points 
is greater than the square of some threshold parameter introduced because of 
the algebraic fall-off of the instanton. However, their methods and results rely 
heavily on the assumptions that the instantons are dilute — an assumption which 
we can see from figs. H and ^| is an unreliable one, particularly at their smaller 
degree of cooling.. 

Our distributions are peaked at around 3/m + at (3 = 2.4, and at 2/m + at 
(3 = 2.5. As in || pl| we find the smaller scale instantons, those with pi < 2a, 



absent after cooling. Thus at smaller lattice spacing, we are able to explore the 
distribution down to smaller instanton sizes without distortion. 



Results presented recently |19| showed a plateau in the distribution of sizes 
calculated on a lattice with twisted boundary conditions; our data, calculated 
using periodic boundary conditions, show no evidence of such a plateau. 



Analagous to the discussion in [0, we can deduce a size distribution of the 



SU (2) vacuum instantons in the dilute instanton gas approximation. In this case 



the distribution is expected to be |2C] 



' ' W // :1 (36) 



VS! dp 



Clearly our distributions do not display this behaviour. It should be noted that in 
this case the dilute instanton gas model gives a size distribution that is infra-red 
divergent, unlike the ultra-violet divergence for 0(3). 



6 Instanton— anti-instanton separations 



As before, in [|TT], we examined the separations of pairs of like and unlike objects 
(I-I, A-A and I-A pairs). There we found that the average closest separation of 
unlike pairs was significantly smaller, with the I-A pairs having a closest sepa- 
ration approximately 70-75% that of the I-I or A-A separations. In section ETT 



above, we showed how the interaction between unlike pairs is attractive and so 
isolated I-A pairs should be found to occur closer on average than I-I or A-A, but 



that I-A pairs closer than a certain minimum distance, given in eq. 20, mutually 
annihilate, so that we should not expect to find I-A pairs at smaller separations 
than when one is located on the boundary of the other with their sizes such that 
-Rmm = Pi = 2p2- While we did not expect to find this relation exactly repro- 
duced by the lattice results, as the instantons are not hyper-spherical and each 
pair will be influenced by every other object in the configuration to a greater or 
lesser degree, we were surprised to find that it is partially reproduced: the closest 
separation of unlike pairs is approximately twice the size of the smaller of the 
pair, but is still significantly greater than the size of the larger. We find that 
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(3 L Method Sep. (R m in) (p c ff (-Rmin)) # disc. # used 



OA ~\ P. Q( rr-\ 

ZA It) >5{X) 


T T 
1-1 


O OA { ~\ 

Z.ZQyLo ) 


Z.i)6{00 ) 


o.ol^Uy ) 


A Q 


1 Q 

zlo 




A A 

A- A 


l.oy^lo ) 


o nn/n^^ 

z.yu^uo ) 






ZlU 




T A 






A 1 7(C\R\ 








U/L 


3.01(26) 










\QW\ 


I-I 


2.02(17) 


2.98(05) 


3.67(08) 


56 


205 




A-A 


1.77(13) 


2.95(05) 


3.68(08) 


57 


204 




I-A 


6.52(09) 


3.18(06) 


4.05(06) 


59 


202 




U/L 


3.44(32) 










2.5 24 \Q(x)\ 


I-I 


1.62(16) 


2.27(04) 


2.84(07) 


44 


206 




A-A 


1.25(11) 


2.21(04) 


2.83(06) 


29 


221 




I-A 


5.92(11) 


2.50(05) 


3.15(06) 


40 


210 




U/L 


4.14(47) 











Table 3: The average closest separation of Like pairs (I-I and A-A) and f/nlike 
pairs (I-A). 'Method' indicates whether the sizes and locations were calculated 
using the action density S(x) or the absolute topological charge density |Q(x)|. 
(peff(-Rmin)) gives the average sizes of the closest objects, averaged over the ensem- 
ble of 250 configurations, with the left figure giving the average size of the smaller 
object. The last two columns indicate how many configurations were discounted 
from the calculations (for example a configuration with only one instanton was 
discounted from the I-I calculations) and how many from the ensemble were used. 
All sizes in units of £. 



there is no appreciable difference in the ratio of size in like or unlike pairs, with 
the larger in the pair being 25-30% as big again as the smaller. Our lattice data 
is shown in table ^[ 

Recently [^l[ it has been proposed that the instantons and anti-instantons 
exist in isolated regions of the lattice, so that 'bubbles' of like objects exist. This 
could account for the observed behaviour in the closest separations shown in 
table [| As a further probe of this we look at the distribution of separations, 
dN/VdR, for both like and unlike objects. Our results are given in figs. |8|-[T0 



The double peak structure of the like separation distribution is evidence that 
objects are localised within groups of like objects. 

This is a striking result and we should consider whether it could be a con- 
sequence of our instanton finding algorithm. If an isolated instanton (or anti- 
instanton) were to have a sufficiently rough peak structure with several subsidiary 
peaks, then this would end up counted more than once. The consistency between 
the result obtained by analysing the S(x) and the |Q(sc)| distributions does sug- 
gest that this is not an explanation of the enhanced signal of like species at small 
separation. 
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Defining a centre-of-mass x c for the (anti-)instantons we looked at the pro- 
portion of objects located within certain radial distances of x c . The results are 
shown in fig. O. The shape of the curve, coupled with the data in figs 



implies that the (anti-)instantons are clustered in many small bubbles, as one 
large cluster would give a very sharp rise in the P(R) curve and a single, broader 
peak in dN/VdR. Clearly this is an area in which further study is warranted. 



7 Conclusions 

Cooling is a local smoothing procedure which enables large scale excitations to 
be studied by reducing the quantum fluctuations. In order to obtain consistent 
results at different lattice spacings, it is essential to calibrate the cooling pro- 
cedure. We explored this thoroughly in 0(3) where we used as a criterion the 
stability under cooling of an isolated instanton of a given physical size. In this 
work, we found that for £77(2) (with the Wilson action), the isolated instanton 
configurations were extremely stable under cooling. Thus, as a calibration, we 
chose instead to create an isolated instanton solution, add the quantum fluctua- 
tions by simulation, and then select the cooling so that it removed these quantum 
fluctuations. Obviously, it would be possible to cool for some arbitrary amount, 
as in the majority of the literature, but we feel there should be some attempt 
made to justify the amount of cooling performed, and this calibration is essential 
if the cooling is performed, and the results compared, at different lattice spacings. 

Our analysis shows that our resulting cooling prescription should preserve 
isolated instantons of size p > a. Of course, cooling could have a more pronounced 
effect in the realistic case of dense instantons: for instance by removing instanton- 
anti-instanton pairs which are close together or overlapping. This is hard to 
calibrate and must be remembered as a possible source of systematic error. 

Compared to previous explorations of topological charge distributions in quen- 
ched £77(2) and SU(3) [[| [Tj|, we have been able to study the topological charge 
distributions in our cooled configurations with large statistics, varying lattice 
spacing and large physical volume. We have made a comprehensive study of the 
size p of the peaks in the topological charge distribution. This size distribution 
increases as p 7 ^ 3 for a dilute instanton gas. We find a similar rapid increase for 
small p but then a maximum for p « 0.3 fm and a rapid decrease at larger p. 
Previous work || gave evidence for a peak in the size distribution on cooled 
configurations, whereas in citechu there is indication of a peak only in their more 
severely cooled data. We have explored this distribution in more detail and at 
different lattice spacings. The absense of important contributions from small 
objects agrees with earlier work. In particular, this explains why cooling gives a 
reasonable estimate |2B| for the topological susceptibility - since the small objects 



which would be modified by cooling are not present anyway. Our conclusions are 
also compatible with phenomenological models for the QCD vacuum WM which 
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require instantons of size p ~ |fm with a density n ~ lfm~4. 

Our results show some evidence that the peak in the size distribution is at 
smaller p for the smaller lattice spacing {(3 = 2.5). We also find a slight difference 
in the instanton density, namely n(/3 = 2.4) = 0.71 fm~4 and n((3 = 2.5) = 0.96 
fm~4. This residual dependence on (3 might be caused by our cooling calibration 
being not quite optimum and so corresponding to slightly stronger cooling at 
(3 = 2.4. Another possible explanation is that lattice artefacts in the instanton 
solution are important for p < 2a and so are more significant at the larger lattice 
spacing. This uncertainty can be resolved by varying the cooling amount and by 
reducing the lattice spacing even further. 

We have studied the relative disposition of instantons (I) and anti-instantons 
(A). We obtain the expected result that there is a minimum IA separation - 
since they will tend to annihilate if too close. We find, however, a surprising 
result: that there is a signal from II and AA pairs at close separation. This could 
arise if instantons were to exist in localised groups, well-separated from localised 
groups of anti-instantons. This clearly needs further work to establish the details 
of such a mechanism. 
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Figures 
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Figure 1: N (S/ Si) for two example SU (2) configurations as they are cooled with 
a = 2. 
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Figure 2: Contour plots for fig. §. S(x, y, 14, 9) is the left plot and Q(x, y, 14, 9) 
is the right plot. 
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Figure 3: The action and topological charge densities for the (x, y, 14, 9) plane 
from a sample configuration, generated at (3 = 2.4 on a 16 4 lattice. 
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Figure 4: Peak height against instanton size p (in units of the glueball mass) 
calculated using S(x) on a 16 4 lattice at f3 — 2.4. 




Figure 5: Peak height against instanton size p (in units of the glueball mass) 
calculated using 
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Figure 6: Size distributions (p in units of the glueball mass) calculated on a 16 4 
lattice at f3 — 2.4. 
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Figure 7: Size distributions (p in units of the glueball mass) calculated on a 24 4 
lattice at f3 — 2.5. 
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Figure 8: The data for the distribution of separations of like and unlike pairs (in 
units of the glueball mass), calculated at f3 — 2.4 on a 16 4 lattice using S(x). 
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Figure 9: The data for the distribution of separations of like and unlike pairs (in 
units of the glueball mass), calculated at f3 — 2.4 on a 16 4 lattice using |Q(x)|. 
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Figure 10: The data for the distribution of separations of like and unlike pairs (in 
units of the glueball mass), calculated at (3 = 2.5 on a 24 4 lattice using 
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Figure 11: The proportion, P(R), of objects found within a distance R (in units 
of the glueball mass) of the centre-of-mass calculated as follows: at f3 — 2.4, 16 4 



using S{x), I (O), A (+), using \Q{x)\, I (□), A (x); and at (3 
|Q(x)|,I(A),A(*). 



2.5, 24 4 using 
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